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0. Introduction 

The Argonne Vi% potential [1] is an updated version of the nonrelativistic Argonne potential that fits both 
np data and pp data, as well as low-energy nn data scattering parameters and deuteron properties. The 
potential was directly fit to the Nijmegen NN scattering database, which contains 1787 pp and 2514 np data 
in the range 0-350 MeV, and has an excellent x 2 per datum of 1.09 . 

The strong interaction part of the potential is projected into an operator format with 18 terms: A 
charge-independent part that has 14 operator components (as in the older Argonne v\a ) 

1, (Ji ■ <jj, Sij, L ■ S, L 2 , L 2 (Ti ■ dj, (L ■ S) 2 

n ■ Tj, (ai ■ aj) (n ■ tj), (n • Tj), L- S (n- Tj), L 2 (n ■ Tj ), L 2 (a l ■ uj) (n ■ Tj), (L ■ S) 2 (n ■ t 3 ) 

And a charge-independence breaking part that has three charge- dependent operators 

Tij , {(Ti • Uj ) Tij , Sij T-ij 

where = 3r zi T Z j — n ■ Tj is the isotensor operator, defined analogous to the Sij operator; and one 
charge-asymmetric operator 

7~zi T z j 

The potential includes also a complete electromagnetic potential, containing Coulomb, Darwin-Foldy, 
vacuum polarization, and magnetic moment terms with finite-size effects. 

Matrix elements of the AW interaction can be specified either in particle-particle (pp) coupling or in 
particle-hole (ph) coupling. Both of the matrix elements completely specify the interaction, and either set 
can be calculated from the other set. We first discuss the symmetries and the angular momentum couplings 
(AMC). 



1. Symmetries of two-body matrix elements 

Generally we use the definition: k = j + 1/2 and j = ^/2j + 1. The matrix element in m-representation can 
be written as 

(m a m b \V\ m c m d ) = V mambiTncmd = (mam^ 1 ^^^^ 1 ) = (m a fh c \V\m d m b ) 1.1 

We use these notations interchangeably. 

The angular momentum coupled matrix element is constructed as 

(Uajb)L \V L \ (j c j d ) L ) = ^2 (j a m a j b m b \LM)(j c m c j d m d \LM)(m a m b \V\m c m d ) 1.2 

ni a mbm c nid 

which is independent of M. It has the symmetry: 

(Uajbh \V L \ Ucjdh) = {-) k « +kMkd ((j b 3a)L \V L \ (j d j c ) L ) 

The inverse relation is 

(m a m b \V\m c m d ) = ^(j a m a j b m b \LM)(j c m c j d m d \LM)((j a j b ) L \V L \(j c j d ) L ) 1.3 

LM 

Instead of using this pp-coupling one can use p/i-coupling to construct angular momentum coupled 
matrix elements. The p/i-coupled matrix clement is defined as 

((jij2)x\V x \(hj 4 ) x ) = Yl R fe2+m2+fc4 ~ m4 0'™'2 - m 2 \X f i){j 3 m 3 j 4 - m^^m^V^ml 1 ) 

1.4 
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These are again independent of /i. 

The matrix elements in pp-coupling are related to the matrix elements in pft,-coupling through [2] 



((^2)A|^ A |(j3j4)A)=E(-) fe2+fe3+i+1 ( 2L + 1 ){^ f, X L }((hh)L\V L \(j 2 j 3 ) L ) 1.5 

This is obtained via general recoupling, using: 

y^^jim 1 j 2 m 2 \LM) {jsmsjjm^LM) = 

M 

= { _y 2 +m 1+H+m3 yw + l)f 31 32 L }<jimij 4 - m i \\ l x){j 3 m 3 j2 - m 2 |A/x) 



A/./. 

={ _y 2 - m2+ji+m3+ L J-(2L + 1){£ j * L x }( Jimi j 3 m 3 |A M )0 4 m 4 . ?2 - m 2 |A M ) 



The inverse relation for the matrix elements is 



{{hk)L\V L \{hh) L ) = (-) L + 1 ^(-) fc2 + fe3 (2A + l){^ & ^}((jlj2)A|^ A |(j3j4)A> 1.6 

Anti-symmetric matrix elements are formed in m-representation as 

(m a m h \V a \m c m d ) = (m a m h \V\m c m d ) - {m a m b \V\m d m c ) 1.7 
In pp angular momentum coupling this becomes 

((jajb)L\V aL \(j c j d ) L ) = ((j a j b ) L \V L \(j c j d ) L ) + (-) L+k < +k «{(j a j b ) L \V L \(j d j c ) L ) 1.8 

These anti-symmetric matrix elements have the symmetry 

((jaj b )L\V aL \(j c j d ) L ) = (-) L+k < +k «((j a j b ) L \V aL \(j d j c ) L ) = (-) L+k « +kb ((3 b Ja)L\V aL \(j cJd ) L ) 1.9 

In p/i-coupling the anti-symmetrization can be obtained by transforming the anti-symmetric matrix elements 
in pp-coupling using (1.5). Thus we obtain 

((jlh)x\V aX \(j 3 J 4 ) X ) = ((jlh)x\V X \(hh)x) + 

!>)A + A'R fe2+fe 3 (2 y + 1) {ii j a A ,}(0,?3)a^ a '|0 2 J 4 )v) L1 ° 

The symmetry of the relations can be significantly improved by introducing the i?m<?-phase. For pp- 
coupling the i?mg-phase is 

((jaj b )L\V aLR \(j c j d ) L ) = (-) k ^((jaJ b )L\V aL \(j cJd ) L ) 1.11 

With this extra phase the matrix elements have the symmetry 

((jaj b ) L \V aLR \(j c j d ) L ) = (-) L ((jaj b ) L \V aLR \(j d j c ) L ) = (-) L ((j b j a ) L \V aLR \(j c j d ) L ) 1.12 

and the anti-symmetric matrix elements can be calculated as 

((jaj b ) L \V aLR \(j c j d ) L ) = ((j a j b ) L \V LR \(j c j d ) L ) + (-) L ((j a j b ) L \V LR \(j d j c ) L ) 1.13 
For matrix elements in p/i-coupling the i?mg-phase is 

((hh)x\V aXR \(j 3 j 4 )x) = (-) kl+k H(jlh)x\V aX \(j 3 j 4 ) X ) 1.14 
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With this phase the anti-symmetric matrix elements in p/i-coupling have the symmetry 
((hh)x\V aXR \(hk)x) = {{hh)x\V aXR \{hh)x) 



v -k h A 

Also, the anti-symmetric matrix elements in p/i-coupling can be constructed as 
((hh)x\V aXR \(j 3 j 4 )x) = ((jih)x\V XR \(hj 4 )x)+ 



-k J, A' 

and the relation between pp-coupling and p/i-coupling is simplified to 

((jij4)L|^ ifl |(j 2 j 3 ) L ) = (-) L+1 X:(2A + l){^ j l}{(jih)x\V aXR \(hk)x) 1.16 

or its reverse relation 

((jlJ 2 )A|^ Afl |(j3j4)A>=E(-) L + 1 ( 2L + 1 ){? f ^}(0U4)L|^ aifl |(j2j3)L) 1.17 

L 

In all the nuclear structure calculations we will use anti-symmetric matrix elements that include the 
i?mg-phase. For simplicity we will drop the superscripts aR throughout. 



2. General considerations and relationships 

We assume the interaction can be written as a scalar product of tensors of rank (k). Then the ph-ph matrix 
element including the -Ring-phase of (— )( fe i+ fc3 ) can be computed as 

< (l,2) x \V H \(3,4)x >=(-)( fel+fe3 ) < (1,2) A M(3,4) A > = 

=( _)(ki+fcs) < (1,2),|[[/«(1)0^(2)]|(3,4) A >= 21 

=S k ,xH kl \ < 1||C/ (A) I|2 > R fc 4 < 4||^ A >||3 > 
A A 

Thus it is necessary to bring the various interactions into this form. 

We will use the Fourier transform to separate the variables fi and r 2 . For this, the general relation 
can be worked out : 

V(r 12 )Yj. M (h2) =- f q 2 dqV J (q) V ^_L(^ 2 0|J0) 

* J £t, J ^ 2.2 

x^ 1 (gr 1 ). 7£2 ( <Z r 2 )( J )(^-^- / )47r[y, 1 (f 1 ) ® Y l2 (f 2 )f M) 
where the form factor of the interaction is 

V J (q) = J V(r 12 )jj(qr 12 )r 2 12 dr 12 2.3 

The various necessary operators are defined as 

ri2 =n-r 2 

S = ^(^i + ^2) 2.4 
L = ^r 12 x (V, - V 2 ) 



Using ([3], eq. 5.1.4) 

and ([3] eq. 5.2.4) 

Using ([3] eq. 5.1.8) we write 



a W Qb (k) = (_)*fc[ a « 0fe «](°) 



=ri2^[YW(ri2) ® (Vi - V 2 ) (1) ] 



2.5 



2.6 



2.7 



(i) 



For terms containing L 2 we need to evaluate the commutator for the tensor components with V 
Vi — V2 which takes the form 

[rf),V«]=25 M 

With this we obtain 



[r <E> V] (1) ® [r <g> V] 



(i) 



1 


1 




1 


1 


il 




Jp 




' 1 


1 






1 





\{Jr 



(3) 



2.8 



lb) 



Due to vanishing a x a, each of j, J r , J p can take up only the values of or 2. For this operator the following 
relation is useful 



[r ® r] ( J) = r 2 ^ [y C 1 ) (f) ® y« (f )] ( J) - r 2 J — ^- < 1010| JO > (f) 



2.9 



Thus we can write 



[r ® V] (1) ® [r ® V 



(i) 



(3) 



=3r 2 ^ J p | 1 1 1 1(1 1 0\J r 0)V4^[y( J ")(f)® [V® V] 

Jr-Jp \ f ^ * 

1 1 11 /4tt 

i 



U) 



-6r{J J [y«(r)®V 



r?l (3) 



2.10 
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3. Central interaction 

Matrix elements for the interaction term V c {r\ 2 ). This term includes the Coulomb interaction. The matrix 
element contributes only for natural parity matrix elements, i.e. A + l\ + l 2 = even. Including the "Ring"- 
phasc of (— ) fe i+ fe 3 wc get 

((12) A |V C >*|(34) A ) = I J q 2 dq F c (q) (^(-) fel <1||Y A ||2) 45 (-P (4||F A ||3)^ 3.1 



where the form factor for the central interaction is 



F c (q) = J V c (r 12 ) j (qr 12 ) r 2 12 dr 12 3.2 

For the evaluation of the matrix element we need the following reduced matrix elements in which we 
include the radial matrix elements involving the Bessel functions: These also contain the Ring-phase. They 
are evaluated for terms that do not contain a using ([3] eq. 7.1.7) 

^£(-) fel <(^i S i)jil|GW||(^ 2 ). 72 >=(-)^ +fel+fc2 ^{ 4 '? , k / AV^<l 1 \\G^\\i 2 > 3.3 
k k W2 Ji L/ti 

We abbreviate this as 

/(fc,l,2)V4^<4||G (fc) p 2 > 3.4 
At this point we need the following reduced matrix elements ([3] eq. 5.4.5): 

V47 < ii\\Y e \\i 2 >= {-)Hdi < hOe 2 0\m > (R^jxiqr)^)) 3.5 

We abbreviate this as 

y(£,l,2)(R 1 (r)\j x (qr)\R 2 (r)) 3.6 
With these definitions we write the reduced matrix elements as 



[n (-)^(l\\Y x \\2) = /(A, 1, 2)y(\,l,2)(Ri (r)\j x (qr)\R 2 (r)) 3.7 



G 



4. Term (LL) 

We write the interaction using eq.(2.8) with j = as 



V = V L2 (r 12 ) [LQ L] = -V3V L2 (r 12 )[LW L (D] W 
= - ^ ^ L2 (r 12 ) [[f 12 ® V12] (1) <8> [ri 2 <8> V12] 



(i) 



(o) 



V3 



f 1 


1 




I— 1 


1 


;} 




Jp 





(«) 



4.1 



+ j J}ri 2 y i2 (r 12 )[r«(r 12 )®V 12 ] 



(0) 



= ^V4^r 2 2 t/ L2 (r 12 )^(1010|J0){| \ j } [r (J) (r 12 ) ® [V12 ® V12] (J) 



(0) 



- ri 2 F L2 (ri2) [y (1) (fi 2 ) ® V12] 



(0) 



As the 9j-coefficicnt requires J r — J p we have set J r = J p = J. J can assume the values of or 2. We 
employ the Fourier transform where we define the form factors of this interaction as 



V 



L2 



= / ri 2 V L2 (r 12 )jj(qr 12 )dr 12 



and similarly 



4.2 



4.3 



V L2 >\q) = j r\ 2 V L2 {r 12 ) 3l {qr l2 )dr l2 
Then, using eq.(2.2) with eq.(4.1) we obtain 

V=J2 §<1010| J0)4{j I )}Y,i-) {tl+l2+J) ' 2 hi2< W*0\ J0 >- I ' q 2 dqV L2J {q) 3t Mri)3 t Mr2) 
x (-Y^nUY^ih) ®Y^{f 2 )] {J) ® [V12 ® V12] (J) 



(0) 



_ <*iM 2 0|10> - / q z dqV^{q)j t Mn)kMr2) 



1(0) 



4.4 

Substituting V12 = Vi — V 2 results in five terms which need to be recoupled in order to be of the form 
F( fc )(l)©G( fe )(2). These are 



Zi =(-) 

1 " 

~T 2 



[Y^\l)®Y^\2)] {J) ® [Vi®Vi 



v4 1 (-0 



0y< <2 '(2) 



(0) 



(o) 



T (£l) (l)® [Vi®Vi] (J) 
Z 2 =(-2)(-Y 2 [[Y (ei) (l) ® T^)(2)] (J) ® [Vi ® V 2 ] (J) 

= 2 E j {t t fc}[[^ (£l) (l)^V 1 ] (fc) 0[T^)(2)®V 2 ] 



(fc) 



<-Y 2 
i 



[y( £i )(i)®y( fe )(2)] (J) ® [v 2 <g>v 2 ] 



V7 1 (-0 



1(0) 



*1 



r (£l) (i)0 r (£2) (2)«) [v 2 <g> v 2 



(^) 



z 4 =(-)" 2 [y (£i) (i)®r^(2)] w ®Vi 
[r (£l) (i)«v 1 ] (fe) 0r( £2 )(2)' 



- l(0) 



1 



4.5 



(_)fer[y(<i)(l)g,y(&)(2)] 

y( £i >(i)© [y^(2)®v ; 



(i) 



i(<i) 



(0) 



For the evaluation of this matrix element we need the additional reduced matrix elements: The matrix 
element involving the operator V: 



4.6 



4ir < £i\\(Y e vY K) \\£2 >=(-f^ K iilh 



v(2£ 2+ 3)(. 2 + d{/ 2 I I ^mrm^-^mr)) 

1 



-i}(o' o 



which we abbreviate as 

^ de/(i,^, «,l,2)<iii(r)|j/(gr) opWI^W) 

i=l,2 

where op(l) = ^ and op(2) = i. 

And the matrix element involving the operator V 2 : 



4.7 



4tt(4|| [V<g>V] (J) 



(A) 



p 2 ) = {-)^ +x+j y\jhi 



A 



:}{ 



1 



1 



J 



£ 2 + 2JU 2 4 + 2 £ 2 + 



,}( 



) 



x (4 + 2)^4 + 1^2 + 2 < Ri(r)\j t (qr)( — 



£2 + In ,_d 

K dr 



1 ± )(l-T)\Mr)> 



A 1 
\ £2 



J A 

h £ 2 i I £2 £2 £2 + 



m £2 £ 2 +i}( 



J \(£\ £ £2 
) . 

x <^<-2 - D< /?i cil.// (y)(^ 



000 

d_ + £2 + >- . 

dr 



./ A | | 1 1 ./ 



£1 £ 2 )l £2 £2 £2 
x£ 2 £ 2 < Ri(r)\j e (qr)( 
1 



\(£i £ £2 
1 i V 

d 



4-l w d 
dr r '^dr 



Is 



\R2(r) > 



£ J A 

£2 £1 £2-2 



2 £2 



i}( 



£1 * 4-2 
00 



x (4 - 2)V£=TV5 < fliWii^^cl + 7)(|: + ^f±)|Jfc(r) > 

4.8 

Here we give the reduced matrix elements separately for J = and for J = 2 since they combine with 
different form factors. For J — the selection rules require £ = A. Thus for J = we can abbreviate (4.8) 
as 

53 dde/ (i,A,l,2) < Ri(r)\j e (qr)op'(i)\R2(r) > 4.9 

i=l,3 



whereas for J = 2 we abbreviate this as 

ddel 2 (i,£, X,l,2) < Ri(r)\j e (qr)oi/(i)\Ra(r) > 

»=1,3 

with op'(l) = £s, op'(2) = if, and op'(3) = 4,. 
Furthermore, we note that 

V°' 00 >{i 11} = ^ 



2* 1 'LI 1 OJ VT2 

and 



The matrix element of the J = piece can be written as 

rR,L2 I / o7\ ^_/_\ki+k 4 ^ [ 9j„frL20t 



< (12) A |V ( ^|(34) A >=(-) fel+fc4 - / <z 2 ^ 20 (g) 

< 1||(^ [V V]<°>)< A >||2 > ^ < 4||F A ||3 > 



4n < i||yA|| 2 > V4E < 4 || (y A [y V ](0)) W|| 3 > 



A A 
This matrix element contributes only in natural parity cases. 

< (12)a|^2 2 (.7=o)K34)a >=v^(-) kl+k4 l J 1 2 dqV L2 °(q) 



E (-)" +A+1 ^ < l||(y / V) A ||2 >< 4||(r £ V) A ||3 > 



<=A-1,.,A+1 ^ 

For the term with J — 1, linear in r, the contributions are from Z4 and Z$: 
<(12) A |^L L2 |(34) A > = 



1 4tt 



( _)(fci+fc4) £ (-)( £ + A+1 )/ 2 £<^0A0|10> - [ q 2 dqV L21 (q) 

e=x-i,x+i 71 •* 

(< i||(r £ v) (A) ||2 >< 4||r A ||3 > + < i!|r A ||2 >< 4||(r £ v) (A) ||3 : 

This term contributes only in natural parity cases. We have to consider £ = A — 1 and £ = A + 1. 
Finally, the contributions from J = 2 can be written 

< (12)a|^L L 2 2 )z=1 |(34) a >= -^^(-) ikl+k * )2 - J <?daV L *>\ q ) 

J2 (-)( £ + A )/ 2 ^<£0A0|20> (< i||[r £ (v®v)( 2 )] (A) ||2><4||y A ||3> 

£=A-2,A,A+2 

1||F A ||2 >< 4|| [Y e (\7 ® V) (2) ] (A) ||3 >) 



+ < 

Again, this term contributes only in natural parity cases with I = A — 2, £ = A, and £ = A + 2. 
The remaining contribution is 

< (12)a|^L L 2 2 )z=2 |(34) a > = 

E E ^ 9 (4,4,A)<i||(F^V)( A )||2><4||(y^v)( A )||3> 

£ a =A±l«(,=A±l 



with 



d lsq (£ a ,e b ,\) = J±(-)V' +t »v 2 iJb < 4040120 > [ l L 1 1} 



4.18 



This term contributes in both, natural and unnatural parity cases. We have to consider the three cases with 
£ a = lb = A — 1, A, and A + 1 as well as the two cases: £ a = X — l,£f, = X + 1 and £ a = X + l,£b = X — 1. 
With our abbreviations the reduced matrix elements can be written as 



47T 



(_)fe 1 ^(l|| ( ^V)A|| 2) = £ f(X,l,2)del(i,£,X,l,2)(Ri(r)\j e (qr) op(z)\R 2 (r)) 4.19 

A i=l,2 

<l||[y A (V®V)(°)] (A) ||2>= £ /(A, 1, 2)ddel (i, A, 1, 2)(i?i(r)|^(qr) qp'(i)|ik(r)) 4.20 

i=l,3 



and 



(_)*! <1||[^( V ®V)( 2 )] (A) ||2>= ^ /(A, 1, 2)ddel 2 (i, £, A, 1, 2)(i?i(r)|^(<jT) op'(z)|i? 2 (r)) 4.21 

i=l,3 



5. Terms aa and (LL)ctct 

Such terms can be obtained from the previous terms such as the central interaction or the L 2 -interaction 
by adding the <J\(7 2 term. This results in 

V = V L2a (r 12 )[LQL][a 1 Q<J 2 ] 5.1 
We can obtain the separated interaction immediately from the separated form for the L 2 term via 



[F^(l) G< J >(2)] [ft ff 2 ] = J2(-) (k+J+1) 

k 



[^ (J) (l)®^i] (fc) [G (J) (2) 



CT2 



|(fe) 



5.2 



Thus, we can immediately take the results of the previous sections and write the matrix element for 
the cr-interaction as 



<(12)a|V ( ^ ) |(34) a >= 



47T 



2A 



_l_\ki+ki'^_ f 
1 7T J 



q A dqV°(q) 



(_)(m+l) < i||(yV) A ||2 >< 4||(yV) A ||3 > 

e=\-i,.,\+i 



5.3 



Similarly we write the contributions to the matrix elements of the L 2 er-intcraction for J = as 

< (12) A |^|(34) A >= * ^-A-f^ - f q 2 dqV L ^(q) £ R^ 1 ) 

V12 2A + 1 n J ,-a-i,.,a+i 



/T2 2A + 

( < 1||[(F £ [V0 V] (0) ) W a] A ||2 ><4||(r £ a) A ||3 > 



5.4 



+ < l||(rV) A ||2 >< 4|| [(Y e [V v] (0) )V 
We include three cases: £ = X — 1, £ = X, and £ — X + 1. 

< (12) a |^=o)I(34)a >=^^T(-) fel+fe4 ^ / ^"""fo) 



13 > 



E E (-) 

K £=B-1,K+1 



(£ + A) 



<H 



(rv) fi (T 



12x41 



5.5 



13 > 
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In this term £ can take the values of A — 2, A — 1, A, A + 1, and A + 2. 

For the term with J = 1, linear in r the contributions are from Z4 and Z§: 



< (12)a|^ 2 ) CT |(34)a >= 



E E (-) (£+K+1)/2 (-) (K+A+1) ^^ T (-) (fcl+fc4) i<^0|10>^ f q 2 dqV L2 ^(q) 

K e= K -i,n+i V3 + 7t J 

(< i|| [(r £ v)< K V] (A) ||2 >< 4||(rv)( A )||3 > + < i||(yv)W||2 >< 4|| [(y £ v)< K V] (A) ||3 >) 

Here k can be A — 1, A, and A + 1. For each case £ can be k + 1 and k—1. 
Finally, the contributions from J = 2 can be written 

< (12)a|^ 2 2 ; z=1 |(34) a >=E(-) (K+A) J20 2A 4 +l ( - )(fcl+fc4 ^ / ^^(l) 
E (-)( £+K )/ 2 )| < £0k0|20 > ( < 1|| [[Y e (V ® V) (2) ] (k) <t] (A) ||2 >< 4||(F k <t)W||3 > 

+ < i||(yv)( A )||2 >< 4|| [[y £ (v ® v)^] (k) <j] (a) ||3 >) 

Again, k can take the values A — 1, A, and A + 1. For each case £ can be k, k — 2, and k + 2. 



<(i2)a|^l l 2 2 ; z=2 i(34) a >-E(-) (k 



+A). 



1 4tt 
7 T2l2A + 



Ji J2 A . 



we abbreviate this as 



Y J 9(W,2)V^<£ 1 \\G^\\£ 2 



5.6 



5.7 



5.8 



(_)(fci+fc*)2 /" ^^,2^ 
1 7T 7 

+ E E ^(4,4,K)<i||[(^v)( K V] (A) ||2><4||[(y^v)( K V] (A) ||3>) 

For the reduced one-body matrix elements containing a we use ([3] eq. 7.1.5) 

( l\ £ 2 £) 

R fel <(£i S i)ji\\(G^a)W\\(£ 2S2 )j 2 >=(-) k \ n hV6{ 1/2 1/2 1 Iv^F < > 5.9 



> 



5.10 



With our abbreviations the reduced matrix elements can be written as 



4ir 



A 



(-)*^<1||(FV) A ||2) =g(£, A, 1,2)2/(^,1,2) (i^r)!^)!^)) 
A 

(-) fcl (l||[[r £ ® V] K ®a] A ||2) = E ff(K,A,l,2)dei(i,^«,l,2)<i2i(r)|j/(gr) op(t)|ifc(r)) 



5.11 



5.12 



i=l,2 



47r (-) fel (l||[[r K (V®V)(°)] (K) ( 7] (A) ||2> = E g(n,\,l,2)ddelo(t,n,l,2)(R 1 (r)\ Je (qr)op'(t)\R 2 (r)) 5.13 

i=l,3 



A 
and 



i7r (-) fe ^l||[[^(V®V)( 2 )] (K) <7] (A) ||2) = E fl(«,A,l,2)ddei2(*,^«,l,2)<fl 1 (r)|j / (gr) p^WI^W) 



i=l,3 



5.14 
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6. Spin-orbit (LS) 



V =V LS (r 12 ) [L« S&] = -V3V LS (r 12 ) [L« S«] 



(o) 



= - r 12 V LS (r 12 )V2^[[Y^(f 12 ) (Vi - V 2 ) (1) ] ( ) 
Wc introduce the form-factor of the LS-intcraction: 



1(0) 



V LS (q) = / r? 2 ^ s (ri 2 )ji(«ri2)dri2 



„3 t/LS 



6.1 



6.2 



and write using eq.(2.2) with J = 1 



q 2 dqV Lb (q)j ei (qr 1 )j e2 {qr 2 ) 



X l X l 2 



<^i0^ 2 0|10> {iPi-i) 



6.3 



x Ait 



'[Y^(h) r^)(f 2 )] (1) (Vi - V 2 ) (1) ] (1) S« 
Or, combining factors 



(0) 



x hi 2 < ^i0^ 2 o|io > 



6.4 



X 47T 



[y^^(fi) F^(f 2 )] (1) (Vi - V 2 )«j (1) (ax + a 2 )« 



(o) 



Inserting S from above, results in four terms as represented in the last line of the previous equation. They 
must be recoupled and brought into the form F^ k \l) G^(2). These terms are: 



Zx = 



[[y^(l)®y £2 (2)] (1) 0V< 1) ] (1) 0ai 1) 



-i(o) 



=-e*(i '; 1} 



y^(2)0 [y £i (i)®V] 



i(fc) 



y f2 (2)0 Vi] (fc) 0ai 

(0) 



(0) 



^2 =(-) 



[y £i (i)«y" 2 (2)] (1) «v (1) 



(i) 



(40 



(1) 



(0) 



[y £l (l)0F £2 (2)] (1) 0V< 1) 

= H £2 V3E{t \ ^[^(^^'^[^(i)®^)] 



6.5 



(fc) 



^4 =(-) 



[y'>(i)®y ,! (2)] (1) ®v (1) 
fc r h h l 



(i) 



(0) 



= H^E^f{T 1 fc} ^(l)0[[^(2)^V 2 ] (fc) 0a 2 
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by introducing 

d so (£ 1 ,£ 2 ,X) = ^£ 1 X<£ 1 0£ 2 0\10>(-)^ +e2+1 ^ 2 {^ ^ \} 
we write the matrix element as 

<(l,2) A |^ s |(3,4) A >=-^-(-)( fe 1 +^ £ -f q 2 d q V LS (q) 

2A + 1 e^ii nJ 

'Y, A, k) < 1|| [[Y e ® V] (fc) ® <r] (A) ||2 >< 4||r A ||3 



dso(£, A, k) < 4|| ® V] W ® o-J ||3 >< i||r A ||2 > 

k 

-d so (£,x,\) < i||(yV)( A >||2><4||(y A v)( A >||3> 
-d so {£, a, a) < 4||(yV)( A >||3 >< i||(r A v)( A )||2 > 

-d ao (*,A,A) < l||(r A ( T) (A) ||2><4||(r £ V)( A )||3> 
-d so (£,X,X) < 4||(FV)( A >||3 >< l\\(Y e V)^\\2 > 



6.6 



6.7 



All the reduced one body matrix elements have been defined before. 

7. Tensor interaction 

We write the tensor operator Si 2 as 

S12 = 3(0^12X0^12) - 7.1 
Recoupling this, S12 can be written as coupled tensor operator 

(II, 

J\[ai® a 2 \ y ' ® [fi2 ® r i2 J v ' 

J 



S12 = 3 Y j[[<ri ® (J) ® [fw ® ri 2 ] (J) ] +V3[a 1 (g> a 2 ] (0) 7.2 



with only J = and J = 2 contributing. In this sum the last term is cancelled against the J = term in 
the sum. The remaining term is 



S 12 = 3 < 1010|20 > W 2 {r 12 ) ® [<ri ® <r 2 ] 



I (2)1 (0) 



7.3 



Separating fi 2 and evaluating the Clebsch-Gordan coefficient (3<1010|20>= \/6) finally gives 
< (1, 2) A | V^|(3, 3)a >= _i^_ ( _ )(fel+fc4 ) £ £ 1 / 

la h 74 

y/6ij b < 4040|20 > ^ b 2 j < i||(y^ <7 )(A)|| 2>< 4 || (r ^ (T )(A)|| 3 > 

where the form factor is 

^ T (?) - y ^ T (ri 2 ). ? 2(gri 2 )r2 2( iri 2 7.5 



13 



8. Term (LS) 2 

This term can be brought into a similar form as the some of the interactions already discussed. We write 



V = V LS2 (r 12 )[L Q ,S][L Q S\ — V LS2 (r 12 )J2 j[[L (1) ® W) ® ® 



(o) 



8.1 



Here j can assume the values 0, 1 and 2. While the terms with j = 0, 1 can immediately be deduced from 
the previous cases, the term with j = 2 needs to be evaluated new. 
For j — we write 



[£®L] (0) ® [S® S^ (0) = J[L0L] [SQS\ = i[L0L](3 + CTi ■ a 2 ) 

o 



8.2 



Thus this part can be added immediately to the previously evaluated interactions by using V L2 - e if 
+ \V LS2 and V L2 ° 
For j = 1 we write 



yL2 + lyLS2 and yL2a,eff = yL2 a + lyLS2_ 



8.3 



This contribution is incorporated by substituting V LS ' e ^ = V LS — jV LS2 . 
Here it remains to evaluate only the matrix elements of 



V = V LS2 {r 12 )Vb\[L^ ® L«] (2) ® [SW ® S«] (2) 



(o) 



8.4 



Using eq.(2.4) we write 



[SW ® S«] (2) = i ( [* ® (2) + 2 ® <? 2 ] (2) + [a 2 ® <? 2 ] (2) ) 

(2) 

Noting that one-body matrix elements vanish for the operator [<ji <S> ci] we write the remaining terms as 



V5 



V = V LS2 (r 12 )^- ® (2) ® [a« ® (2) 



(0) 



8.5 



/9 s ! 

We assume that similar to eq.(4.4) the operator [L ® L] can be written as a tensor product of variables 
fi and f 2 as [C/( Kl )(l) ® (2)] (2) Thus we can write 



[L®i] (2) ® [ft ® ft] (2) ] (0) - [[[/( K1 )(1)®VF^)(2)] ® [< 



T ® (T 2 1 (2) 



(0) 



^EH^Ii * fc}[[^ (Kl) (l)^-i] (fc) 0[^ (K2) (2)®a 2 ] (fe) 



(0) 



8.6 



We now have to factorize the operator [L ® L] using eq.(1.10) with j = 2. Again, as they combine with 
different form factors, we list the three cases separately namely: (a) J r = 0, J p = 2, (b) J r = 2, J p = 0, and 
(c) J r = 2, jp = 2. Numerically evaluating the 9j-symbol and CG-coefficient we find 



^ [L ® L\ (2) = V^(^r 2 [yW (f) ® [V ® V] (2) 

r (2) (f)® [v® v] (2) 



+ 



VE 

L2 ? 
14 



12 



(2) _ \/2 2 

24 r 
(2) 1 /5 



r( 2 )(f)® [v®v] (0) 

j/fr[y (1) (r)®V] (2) ) 



(2) 



5.7 
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In the following we list all the eleven terms that contribute to [L <g> L] in their recoupled forms: 

z x = [ [y 1 - (i) ® y £2 (2)] (2) ® [V! ® V!] (0) ] (2) 

= "[y £l ® [Vi®Vi] (0) ] (/l) ®r /2 (2)" 
z 2 = - 2 [[y £i (i) ® y fe (2)] (2) ® [Vi ® v 2 ] ' 



(2) 



1(0) 



(2) 



2 

71 



jH^i,^ £ J}[[y^(i)Vi] (Ll) ®[y^(2)v a ] 



(^2) 



(2) 



(0) 



y £i (i)® [y £2 (2)® [v 2 ® v 2 ] 
z 4 =[[r £l (i)®^ 2 (2)] ( % [Vi®Vi] 



(<2) 

(2) 



(2) 
(2) 



® y £2 (2) 



1(2) 



Zr, = 



[r £i (i)® [V!®V!] (2) ' 
-2[[y^(i)®y^(2)] (0) ® [V!®v 2 ] (2) 

= %, (2 EH (il+4) x{l 2 Li ^[[^Vi^M^V,] 
: ^(-)(^)^ 



(2) 

2 



(2) 



?i\/5 



y^(i)® [y £2 (2)® [v 2 ®v 2 



(2) 



(fc) 



(2) 



1(2) 



(2) 



Z 7 =[[Y^{l)®Y k {2)} (2) ® [ViOVi] 

= V5Y / (-f 1+k) k{z I J}[[y £l (i)® [Vi®v 1 ] (2) ] (fe) ®y £2 (2) 

fc 2 1 

-2[[y £i (i)®y £2 (2)] (2) ® [Vi®V 2 ] 



(2) 



z* = 



I (2) 



l(2) 



^2 2 



= -10 ^ L^J 1 1 2 [ [y £i (i)Vi 



(il) 



L x L 2 2 



[y f2 (2)v 2 ] 



(■£-2 



(2) 



fc ! 2 '■ 
?io=[[y /l 



(2Mfc) 



(2) 



y £2 l (1) ®Vi 



(2) 

1 2 1 



5»<* + *>V3*{£ £ ^}[[y^(i)v 1 ] (fe) ®y £2 (2) 



(2) 



(A') 



(2) 



Combining this with with eq.(8.7), we can write the matrix element contributions for the various form 
factors (J = 0, J = 1, and J = 2) as 

\/5 47T 



< (12) A |V« L L f |(34) A > 



... 



24 2A + 

( < i||[(y £ ® [v®v] (2 V K) ^] (A) l|2 ><4||(yV)( A >||3 > 
i||(yV) (A) l|2 ><4||[(y £ ® [v®v] (2 V K) <?] (A) l|3 >) 
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8.9 



+ < 



} 



8.10 



2 - J q 2 dqV LS20 (q)x < 1|| [(^V)^a] (A) ||2 X 4|| [(F^V)^a] (A) ||3 >) 
The contributions from terms Zio, and Zu, the J = 1 terms are 

5 4v H (ii+y E^<«iio>{; i Da 1 1} 



<(12) A |t/«L i 1 f|(34) A >= 



12 2A+1 



1||(FV)( A )||2 ><4||[(r £ »V)^a] (A) ||3 >) 



> 



+ < 

Finally, the J = 2 terms are from Z\, Z 2 , Z3, Z7, Zg, and Zg. 

47T 



8.11 



< (12) A |^L L 2 f |(34) A > 



2A + 1 nJ £7% 

JV^W 1 1 2 U 2 2 2 lr-}("+*) 

12^ 24 U h A/Ui ^ k/ 1 j 

( < 1||[[F £2 <g> [V® V] (2) ] (K) ^] (A) ||2 >< 4||(r £l a)( A )||3 > 
+ < l||(y £l CT) (A) ||2 >< 4|| \\Y l2 ® [v® v 

(A) 



+fe)/2 



8.12 



l(2)-|(«). 



(A), 



||3> 



< (12)A|^ 2 f |(34) A > 



I a}( < [V® V] (0) ct] W ||2 ><4||(r £l o ; ) (A) ||3 

+ < i||(r £i ( ?)( A )||2><4||[[y £2 ® [V®V] (0) a 

_(_)(fci+fc*) 2 /" q 2 dq v^( q ) y ij 2 < £ l0 ^ 2 0|20 > (-)(*+* 

1 •/ ri_ 



13 > 



47T 

2A + 



2)/2 



+ Vl4-(-) 



4 £2 2 • 
1 1 2 
.4 4 2. 



< 1|| V)^<r] (A) ||2 >< 4|| [(Y^V)^a] W ||3 > 



8.13 



9. The isospin component 

The isospin dependence of the matrix elements is best worked out using the (jap) coupled matrix elements. 
In that case all we have to do is to work out the expectation value of the isospin operator. We do this here 
for the four isospin operators that appear in the Argonne vis potential of Wiringa et al [1] and the four 
possible matrix elements (1, tit 2 , T 12 = 3t z it z2 — tit 2 , and t z \ + t z2 ): 



< (pp)lll(pp) > = 


+ 1 


) 


< (pp)\t z i +T z2 \(pp) > 


= +2 


< (pn) 1 \{pn) > = 


+ 1 


) 


< (pn)\T zl + T z2 \(pn) > 





< (pn) 1 \(np) > = 





) 


< (pn)\r z i + r z2 \(np) > 





< (nn)\\\(nn) > = 


+ 1 


) 


< (nn)\r z i + r z2 \(nn) > 


= -2 


< (pp)\nT 2 \(pp) > 




+ 1 


, < (pp)\Ti 2 \( PP ) > 


= +2 


< (pn)|fif 2 |(pn) > 




-1 


< (pn)\T 12 \(pn) > 


= -2 


< (pn)|fif 2 |(np) > 




+2 


< (pn)\T 12 \(np) > 


= -2 


< (nn)|fir 2 |(nn) > 




+ 1 


< (nn)\T 12 \(nn) > 


= +2 
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Thus we calculate the matrix elements 
for < (pp)\V\(pp) > 

< (12)|V|(34) >=< (12)|V + V TT + 2V tT + 2V TZ + V Coul \(34) > 

for < (pn)\V\(pn) > 

< (12)|V|(34) >=< (12)|V - V TT - 2V tT \{34) > 

for < (pn)\V\(np) > 

< (12)|V|(34) >=< (12)|2y TT - 2F tT |(34) > 

and for < (nn)\V\{nn) > 

< (12)|V|(34) >=< (12)|V + V TT + 2V tT - 2V TZ \(M) > 
10. CM-correction terms 

From the center-of-mass (CM) corrections we have two additional terms that must be included in the 
two-body matrix elements of the interaction. These are: 

Pi - Pi _ | Ti 2 Vi ■ V 2 
mA mA 

and 

n ■ r 2 



mA 

These can be computed using eq.(2.1). These matrix elements contribute only for A = 1 and can be 
computed with eq.(3.3) resulting in 

- JL < 12^134 >= < HIVII2 > (-)^-L < 4||V||3 > 



and 

-L < m r T , AU >= _L ( _ )fa _L < m2 > ,_ )fa _L < 4 „ 3 > 

Here 

(-) fe ^ < 1||V||2 >= ^ < h ^2^110 >< Ri(r)\^ + \{h -h + l)\\Mr) > 



and 



,fci 1 .iii„-iio. iiis i, -i, 



< l||r1|2 >= < — 110 >< flxWIrl^W > 



11. The radial integrals 

We assume that each wave function is expanded in harmonic oscillator functions corresponding to the 
oscillator length parameter ok. 

R i , l (r) = Y, A iH n ,i(r) 10.1 

n 

As such each of the radial wave functions can be written as a polynomial multiplied by a gaussian of 
argument — (V2r/ok) 2 . We now introduce the variable x — \f2rjok. With this, each radial wave function 
can be written in terms of x as 

Ri,i{x) =R i ,i{x)e- x2 l i 10.2 
17 



where Rij is a ploynomial in x. Generally, we define the barred functions as those without the exponential. 
In terms of x the wave functions are normalized such that 

R 2 l {x)x 2 dx = (^-^j 10 - 3 

In order to carry out the Fourier transform we use the fact that the Fourier transforms of the harmonic 
oscillator functions are again harmonic oscillator functions in fc-space. Thus we expand the product of two 
radial functions as 

R lM {x)R oh (x) =Y,Bti j H ntL {x) 10.4 

n 

Using Gaussian integration, the expansion coefficients can be written as 

=J2H n , L {xk)Ri,i i (xk)R j ,i j (x k )xlwk 10.5 
k 

with this, we write the integral 

/ Ri, h (r)j L (qr)R jtlj (r)r 2 dr = J]B«ff ni£ (g) 10.6 

where q — \f2ok q. The radial integrals can now be computed according to 

q 2 dq V(q)(R iM \j L (qr)\R Uj )(R Sil3 \j K (qr)\R tJt ) = 

= E E B n B ™ \ I ^ V{q)H n>L {q)H mtK (q) 1 . 7 

n m 

=■ B i,j B s t J L K 

' / j / j n rn n , m 

n m 

Inserting the expression for the expansion coefficients B we find 



q 2 dq V(g)<^|jL(^)|^><iia^b>(gr)|^t> 



■3,h i 

2 _ _ 

= - E E Ri ' 1 * ( Xk ) R hk ( X k)x 2 k W k R s J s (x q )R t ,l t (x q )x 2 q W q ^ 7 n,m H n,L {x k )H m , K (x q ) 
k q m.n 

k q 

10.8 

With this notation all terms can be computed as 

< (12) A |V tot |(34) A >=Y J Ri(x l )R2(x l )Sl 3 R 3 {x j )R i (x J ) 



+Ri { Xl )R' 2 ( Xi )S?j R 3 {xj)Ri { Xj ) 
+R 1 {x^ f > <~ w3 5 '' 
+R 1 (x 
+R 1 (x 
+Ri(x 



)R 2 (x i )Sl j R' 3 (x j )R 4 (x j ) 10.9 
)R' 2 , (x i )StjR 3 (x j )R 4 (x j ) 
)R 2 (x i )Sl j RUx j )R4(x j ) 
)R' 2 (x t )Sl ] R' 3 (x 3 )R 4 (x J ) 



18 



Correspondingly, the exchange contribution can be written as 

< (l2) x \V tot \(M) x >=Y,Mxi)Mxi)El j R2(x j )R 4 (x j ) 
+R 1 (x i )R' 3 (x i )El j R 2 {x j )R 4 {x j ) 

+R 1 (x i )R 3 (x i )El j R' 2 (x j )R i (x j ) i .lO 
+R 1 (x i )R' 3 \x i )El j R 2 (x j )R4(x j ) 
+R 1 (x i )R 3 (x i )El j R^{x j )R 4 (x j ) 
+R 1 (x i )R' 3 (x i )El j R^(x j )R 4 (x j ) 

As the marices S and E only depend on the angular momenta involved but not on the quantum numbers 
rii,ri2,ri3, and n 4l the strategy is to compute the matrices first and then compute all the matrix elements 
with the same angular momentum combinations. This way a large number of matrix elements can be 
computed in which orbits differ only by the quantum number n. 

Note. This strategy is essential in computing the G-matrix interaction in calculating the term: 

J2 J2 Z P3P,Mh 2 < (P3P±)k\V\(p iP 2)k > 

K P3,P4 
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